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THE  PROBABILITY  OF  HITTING  A  POLYGONAL  TARGET 


1.  INTRODUCTION 

Military  targets  are  generally  represented  by  simple  geometric 
shapes  in  performance  models  which  require  the  conputation  of  hit  probabil¬ 
ity.  Circles  and  especially  rectangles  whose  sides  are  parallel  to  the 
coordinate  axes  of  the  weapon  delivery  error  distribution  are  the  most 
commonly  used  shapes.  A  much  more  realistic  representaion  of  target  shape 
would  frequently  be  possible  through  the  use  of  a  polygon.  However,  with 
such  a  representation  the  calculation  of  hit  probability  becomes  much  more 
difficult.  The  purpose  of  this  report  is  to  present  a  simple  algorithm  for 
the  conputation  of  such  probabilities. 

2.  DISCUSSION 

The  following  paragraphs  state  the  problem,  the  assumptions  made, 
and  develop  the  conputational  algorithm. 

2.1  Assumptions. 

The  target  representation  will  be  assumed  to  be  a  general  n-sided 
planar  polygon.  The  weapon  delivery  error  will  be  assumed  to  be  bivariate 
in  the  plane  of  the  target.  The  general  algorithm  to  be  presented  does  not 
require  that  the  delivery  errors  be  normally  distributed,  but  the  FORTRAN 
subroutine  and  numerical  examples  included  will  be  based  on  that  further 
assunption. 

2.2  Problem  Statement. 

The  problem  is  to  derive  an  algorithm  which  conputes  P^,  the 
probability  that  a  random  shot  from  the  assumed  bivariate  delivery  error 
distribution  hits  the  polygonal  target. 

2.3  Method  of  Solution. 

The  basic  building  block  for  this  algorithm  is  the  integral  of 
the  probability  density  function  of  delivery  errors  over  a  region  consisting 
of  the  entire  target  plane  "below"  a  specified  line  segment.  Specifically, 
given  a  line  segment  connecting  point  i  whose  coordinates  are  (xj,y,)  with 
point  j  whose  coordinates  are  (xj.yj),  we  will  need  to  integrate  the  delivery 
error  distribution  over  the  region  R  defined  by 

f  x  between  Xj  and  Xj 

i)  . 


R: 


y  less  than  y^  + 


(1) 


Such  a  region  is  shown  as  the  shaded  area  in  Figure  1.  If  f(x,y)  denotes 
the  probability  density  function  for  the  bivariate  distribution  of  weapon 
delivery  error,  this  integral  is  given  by 


(2) 


This  integral  will  be  positive  when  Xj>Xi,  and  in  that  case  will  be  the 
probability  of  hitting  the  region  R.  However,  in  the  present  application, 
this  inequality  will  not  always  hold;  so  we  will  not,  in  general,  refer  to 
this  integral  as  a  probability. 

Applying  the  rule  for  interchanging  the  upper  and  lower  limits  of 
integration,  it  follows  that 


This  relationship  will  be  important  in  simplifying  subsequent  results. 

The  next  step  is  to  determine  how  these  Fij's  can  be  used  to 
determine  the  probability  of  hitting  a  triangular  target.  The  result  for 
the  triangle  will  be  used  in  a  mathematical  induction  development  to  extend 
the  result  to  a  general  n-sided  polygon. 

There  are  two  possible  types  of  triangles  that  must  be  considered. 
In  the  first  of  these  (Figure  2a),  the  "upper  boundary"  consists  of  two 
sides,  and  the  "lower"  boundary  consists  of  one  side.  In  the  second  type 
(Figure  2b),  the  "upper"  boundary  consists  of  one  side  and  the  "lower" 
boundary  consists  of  two  sides.  The  probabilities  of  hitting  these  two 
types  of  triangles  will  be  derived  in  terms  of  the  F  function,  and  it  will 
be  shown  that  the  sane  formal  expression  can  be  used  for  each. 

For  triangular  targets,  the  hit  probability  is  the  probability  of 
hitting  between  the  "upper"  boundary  and  the  "lower"  boundary.  This  can  be 
•ixpressed  as  the  probability  of  hitting  below  the  "upper"  boundary,  minus  the 
,>r  lability  of  hitting  below  the  "lower"  boundary.  In  the  first  type  of 
triangle  this  probability  is: 

PH(TyPe  d)  =  F1 ,2  +  f2,3  -  F1 ,3  *  <4) 

where  Fj^  +  F2t3  is  the  probability  of  hitting  below  the  "upper"  boundary, 
and  Fi  31s  the* probability  of  hitting  below  the  "lower"  boundary,  with  the 
vertices  labeled  consecutively  in  the  clockwise  direction  around  the 
triangle.  Similarly,  for  the  second  type  of  triangle. 
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Ph  (Type  b)  *  Fi  2  ■  ( Fi  ,3  +  F3  2)  =  ^1,2  "  ^1,3  "  ^3  2  (5) 

where  Fjf2  is  the  probability  of  hitting  fieiow  the  “upper*  boundary  and 
Pi, 3  +  P3  2  is  the  probability  of  hitting  below  the  "lower"  boundary. 

Applying  £ne  relation  Fi  j  =  -Fjj  (Equation  3)  to  each  of  the  "lower" 
boundary  terms  in  these  £wo  expressions, 

PH  (Type  a)  =  Fi  2  +  f2.3  +  F3  1  (6) 

and 

Ph  ( Type  b)  =  Fjf2  +  F2,3  +  F3J  .  (7) 

Thus,  the  two  types  of  triangles  can  be  treated  the  same  formally:  the 
sunmation  of  the  F's  proceeds  in  a  cyclic  fashion  around  the  triangle, 
taking  the  sides  in  turn  and  integrating  in  each  case  from  the  first  value 
of  x  to  the  next  value  in  the  cyclic  order.  The  relative  magnitudes  of  the 
x's  themselves  will  automatically  determine  the  proper  sign  for  the  associated 
value  of  F,  so  that  the  F's  can  always  be  added  to  obtain  the  desired 
result.  Therefore,  for  any  triangle  target  with  the  vertices  ordered  con¬ 
secutively  in  a  clockwise  order,  the  hit  probability  is 

pH  (3)  =  Fl ,2  +  f2,3  +  F3 ti  (B) 

This  result  can  be  rewritten  as 

3 

PH  (3)  =  l  Fj  imod3+l  (9) 

i  =1 

where  imod  3  denotes  the  remainder  when  i  is  divided  by  3. 

Actually  desired  is  the  extension  of  this  result  to  a  simple  n- 
sided  polygon,  that  is,  one  whose  sides  never  cross  each  other.  Since  such 
a  polygon  also  has  n  vertices,  the  desired  extension  is: 

n 

pH(n)  =  I  F-jjmodn+1  (10) 

i=l 

This  result  will  be  established  using  mathematical  induction.  In 
particular,  it  iiust  be  shown  that: 

I:  The  result  holds  for  n=3,  the  polygon  with  the  least  possible 
nimber  of  sides,  and 

II:  Assuming  the  result  true  for  a  polygon  of  (k-1)  sides,  show 
that  this  leads  to  its  truth  for  a  polygon  of  k  sides. 

Part  I  has  already  been  demonstrated,  leaving  only  Part  II.  At  this  point 
t*o  additional  propositions  are  required.  These  are 
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•  Any  k-sided  polygon  can  be  partitioned  into  a  triangle  and  a 
polygon  of  (k-1)  sides,  and 

•  Given  two  polygons  that  share  a  common  side,  a  traversing  of  the 
perimeter  of  each  polygon  in  the  clockwise  direction  results  in 
the  common  side  being  traversed  in  one  direction  as  one  polygon 
is  being  traversed,  and  in  the  opposite  direction  as  the  other 
polygon  is  being  traversed. 

In  "Analytic  Function  Theory,  Vol.  1"  by  Einar  Hille,  Blaisdell 
Publishing  Company,  1959,  there  is  a  proof  (page  286)  that  every  simple 
closed  polygon  can  be  triangulated,  that  is,  partitioned  into  triangles. 

This  result  establishes  proposition  (1).  Figure  3  shows  a  nunber  of  poly¬ 
gons,  some  more  odd  shaped  than  others,  with  each  partitioned  into  a  triangle 
and  a  polygon  of  one  less  side.  Figure  4  illustrates  the  idea  of  the  side 
oonrnon  to  a  pair  of  polygons  being  traversed  in  opposite  directions  as  each 
polygon  is  traversed  in  a  clockwise  motion.  A  proof  of  proposition  (2)  has 
not  been  found,  but  its  truth  is  almost  obvious,  and  will  be  assumed. 

The  proof  of  Part  II  of  the  mathematical  induction  proceeds  as 
follows:  Partition  the  k-sided  polygon  into  a  triangle  and  a  polygon  of 
(k-1)  sides.  This  is  possible  by  proposition  (1).  This  is  accomplished  by 
connecting  a  pair  of  vertices  of  the  k-sided  polygon  which  are  separated  by 
a  single  vertex.  Denote  the  vertices  connected  in  this  way  as  vertices  m 
and  im-2. 


By  the  assumption  of  Part  II  of  the  induction,  the  probability  of 
hitting  the  (k-1)  sided  polygon  is 

Pn(k-l)  =  Fl,2+F2,3+**»+Fm-l,m+Fm,m+2+Fm+2,m*3+*‘*+Fk,l  •  (J1) 

By  the  result  of  Part  I  of  the  induction,  the  probability  of 
hitting  the  triangle  is 

fh(3)  =  Fm,m+1  +  Fm+l,m+2  +  Fm+2,m  •  (^) 

Proposition  (2)  insures  that  the  side  common  to  the  (k-1)- sided  polygon  and 
the  triangle  (the  side  connecting  vertices  m  and  m+2)  will  be  traversed  in 
opposite  directions  in  the  development  of  Pn(k-l)  and  Ph(3).  Thus,  when 
these  two  probabil ities  are  added  to  obtain  Pfl(k),  Equation  (3)  insures 
that  Fm,irH-2  an<*  Frm-2,m  cancel  each  other,  establishing  the  desired  result, 
namely.’that 

pH(k)  =  F1 ,2+f2 ,3+* • ,+Fm-l ,m+Fm,m+l+Fm+l ,m+2  + . +Fk,l  •  (13) 

k 

=  I  Fi,imocJc+l 
1=1 

This  completes  the  proof  that  establishes  Equation  (10).  This 
result  depends  only  on  the  coordinates  of  the  vertices  of  the  polygonal 
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target.  It  becomes  extremely  easy  to  apply,  requiring  only  a  means  for 
calculating  the  function  F-ji  defined  in  Equation  (2).  When  presented  with 
a  general  polygonal  target.’the  vertices  are  numbered  consecutively  in  a 
clockwise  rotation  around  the  polygon,  starting  at  an  arbitrarily  chosen 
vertex.  The  values  of  the  F,j  function  for  the  sides,  taken  in  order,  are 
then  summed  to  obtain  the  hit ’probability. 

3.  COMPUTATIONAL  FORM 

The  basis  for  this  hit  probability  formulation  is  the  function 
Fj t j  defined  in  Equation  (2),  and  repeated  here  as  Equation  (14). 

y (*-«,) 


y  —  -CO 

Weapon  delivery  error  is  generally  assumed  to  be  distributed  according  to 
an  uncorrelated  bivariate  normal  distribution.  Under  this  assumption. 


j  1 

f(x,y)  dy  dx 


(14) 


i»J 


X=Xi 


f(x,y)  =  - 

°x°y 


(15) 


where  (A,B)  is  the  mean  impact  point,  ox  and  cy  are  the  standard 
deviations  of  the  delivery  error  in  the  x  and  y  directions  respectively, 
and 


4>(t)  = 


e 


for  all  t,  the  unit  normal  probability  density  function.  Thus 


X=Xi 


(lb) 
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x-A  y-B 

Make  the  changes  of  variable  u  =  _  and  v  = _ ,  and  let 

°x  °y 


m 


ij 

written  as 


yj-Yi  Xi-A  xj-A  yj-B 

=  -  ,  Sj  =  - and  Ij  =  - ;  then  this  integral  can  be 

0* 


xj-Xj  °x  °x 


Fi,j  =  /J4>(u)*(Ti  +  — -nijjtu-SjlJdu 

U  =  S-j  Oy 


where  <t>(Z)  =  /  <j>(t)dt,  for  all  Z,  the  cumulative  normal  probability 


t=-«® 


function. 


In  general,  this  integrat’''-’  must  be  performed  numerically.  How¬ 
ever,  there  are  two  special  cases  wh.re  this  is  not  necessary.  In  the 
first,  identified  by  Xj  =  xj  (or  Sj=Sj),  the  upper  and  lower  limits  of 
integration  are  the  same  and  Fj  j  =  0.  This  is  consistent  with  the  notion 
that  Fj j  is  the  integral  of  the  probability  density  function  "under"  the 
line  segment  connecting  (xj,y<)  and  (xj,yj).  In  this  special  case,  the 
connecting  segment  is  vertical,  and  thus  has  no  neasurable  area  "under"  it. 
In  the  second  special  case,  identified  by  yj=yj  (or  aiij=0),  the  integral 
reduces  to 


Sj 

Fi,j  =  /  4>(u)#(Ti)du  =  4(Ti)0(Sj)-4>(Si)] 

u=Sj  (la) 

The  function  *  is  available  on  most  computers.  In  the  general 
case  requiring  numerical  integration  to  evaluate  Fj  j,  one  further  simpli¬ 
fication  is  useful.  The  function  4>(u),  is  essentially  zero-valued  for 
;uj>4.  Thus,  the  numerical  integral :-on  can  bn  unde  .tore  efficient  by 
adjusting  the  values  of  Sj  and  $;  so  that  no  integration  steps  outside  the 
i  ltnr/jl  (-4<  u  <  4)  are  included. 

A  FORTRAN  subroutine  For  calculating  the  probability  of  hitting  a 
mV/j  nal  target  using  this  'iiethod,  along  with  instructions  as  to  how  to 
it;  .i  is  subroutine  in  a  FORTRAN  projran,  are  iuolidid  in  the  Appendix. 

4.  SAVPi.c  CALCULATIONS 

The  polygonal  target  shown  in  Figure  5  was  used  as  a  basis  for 
checking  this  procedure  and  the  associated  FORTRAN  subroutine.  Hie  proba¬ 
bility  of  hitting  this  polygon  was  easily  hand-calcul ated  by  partitioning 
it  into  rectangles  each  of  whose  sides  are  parallel  to  the  coordinate  axes. 
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Figure  5.  Polygon  (or  Sample  Computation. 


Figure  6.  Polygon  Rotated  30* 


m 


This  result  was  then  compared  to  the  result  from  the  approach  of  this  report. 
As  the  polygon  appears  in  Figure  5,  the  two  special  cases  of  the 
method  apply  -  for  the  vertical  sides  (such  as  the  side  connecting  vertex 
11  and  vertex  12)  the  function  Fjt j=Q  and  for  the  horizontal  sides 
(such  as  the  side  connecting  vertex  10  and  vertex  11)  Equation  (18)  applies. 
The  agreement  between  the  hand  conputation  and  the  subroutine  result  was 
established.  Then  the  polygon  was  rotated  30°  about  the  origin,  as  shorn 
in  Figure  6.  As  long  as  <Jx=<*y  anC*  A*B=0,  the  probability  of  hitting 
the  rotated  polygon  should  beJ the  same  as  that  for  the  unrotated  polygon. 

In  this  case,  however,  the  general  expression  for  Fjj  applies  for  all 
sides  since  none  of  them  are  parallel  to  the  coordinate  axes.  The  resulting 
probability  was  the  sane  as  that  for  the  unrotated  polygon.  Figure  7  shows 
the  probability  of  hitting  this  rotated  polygon  for  A=B=0  and  as  a  function 
of  0x=oy* 
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LITY  OF  HIT 


0  .2  .4  .6  .8  1.0  1.2  1.4 

°ir  °"y 

Figure  7  Probability  of  Hitting  Rotated  Polygon. 
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APPENDIX 


FORTRAN  FUNCTIONS 

This  appendix  presents  the  following  three  FORTRAN  functions 
required  for  the  use  of  this  method  for  computing  the  probability  of  hitting 
a  polygonal  target.  These  functions  are 

1.  PHIT  (called  by  user's  main  program) 

2.  FIJ  (called  by  PHIT) 

3.  FUN  (called  by  FIJ) 

The  user  mist  include  all  three  of  these  in  his  main  program,  but  needs 

only  to  call  PHIT.  The  arguments  In  the  call  statement  are 

N  =  number  of  sides  in  the  polygonal  target 

X  =  array  name  holding  X  coordinates  of  vertices  of  polygon 

Y  =  array  name  holding  Y  coordinates  of  vertices  of  polygon 

XB  *  X  coordinate  of  mean  of  hit  distribution 

YB  =  Y  coordinate  of  mean  of  hit  distribution 

SX  *  standard  deviation  of  hit  distribution  in  x 

SY  »  standard  deviation  of  hit  distribution  in  y 

The  size  of  the  x  and  y  arrays  (>N)  nust  be  established  with  a  DIMENSION 
statement  In  the  calling  program. 

The  actual  numerical  integration  is  acconpl i shed  with  a  Simson 
integration  subroutine  (SIMSON),  which  requires  the  function  FUN  to  calcu¬ 
late  values  of  the  Integrand.  SIMSON  is  resident  on  the  CDC  computer  at 
Aberdeen  Proving  Ground,  but  users  of  other  computers  may  have  to  substitute 
another  integration  subroutine.  In  the  subroutine  used  here,  a  relative 
error  of  0.001  has  been  specified.  This  should  provide  values  of  Fj  j 
accurate  to  the  third  decimal  place.  The  accuracy  of  the  hit  probabilities 
depend  on  the  number  of  these  F's  that  iwst  be  added,  but  in  general  should 
be  good  to  at  least  two  places,  and  probably  more. 
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FUNCTION  PMIT (NjX/V#  Xt#  YB#SX#  ST) 

*** 

THIS  FUNCTION  COMMUTES  THE  PROBABILITY  0*  HITTING  AN 
N-SIDEM  POLYGONAL  TARGET.  THE  COORDINATE  $  nr  THE  VCRTICFS 
OF  THE  POLYGON  ARE  STORED  IN  THE  X  AND  v  ARRAYS  IN 
CONSECUTIVE  CLOCKWISE  ORDER  AROUND  THE  POLYGON.  THE 
DELIVERY  ERROR  DISTRIBUTION  IS  ASSUMED  TO  BE  BIVARIATE 
NORMAL  WITH  MEAN  ( XB* YB )  AND  WITH  STANDARD  DEVIATIONS 
SX  AND  SY  RESPECTIVELY. 

♦  ** 

DIMENSION  X(1)#Y(1) 

1  FORMAT! »  BOTH  SX  AND  Sv  MUST  DC  GPCATER  THAN  ZERO'/ 

*  '  S<  •  SF1J.3,'  SY  -  •  #  F 10. 3 ) 

IF (SX.G^.O.  .AND.  SY.GT.O.)GOTO  5 
"RIMT  l,  SXjSY 
ST 

5  PHIT-O. 

DO  10  I«1»N 
J-M  JDC#N)P1 

"!iiT-''iiiT*FiJ(x<n,Y(:),x(j)#Y(  jj,xn,Yf.,sx,sv) 

10  CONTINUE 
RETURN 
END 


FUNCTION  riJ(XI/YI#XJ/YJ# XB*YD»  SX. SY) 

«** 

THIS  FHMCTIOM  PERFORMS  THE  INTEGRATION  NT  THE  b I VAF lATr 
HIP.  IAL  3CNSITY  FUNCTIOH  •  U»! Dr R •  THT  DIRTCTFO  LINE  S  EOMr NT 
FROM  P  IINT  ( X  It  YI )  TO  THE  "01 NT  <XJ»VJ). 

**♦ 

COM  DM  /ABC/  A#  B#  SI 
EXTERNAL  FUN 
FIJ-T. 

IF  (XI.EQ.XJ)RETURN 

cxi-xbj/sx 

SIP-SI 

si"-anax:<si*#-*.) 

5  IP-AMIH*!  S  I°#<.» ) 

SJ- (XJ-<B)/SX 
SJ^-SJ 

SJP-AMAXl(SJP,-«,.) 

SJR-AMINKSJP#  A.) 

IFCI®.E1.SJ’,)RETHRN 
/«■!  YI-YO)  /SY 

D«SX*!VJ'-YI)/!$Y*(XJ-XI)) 

inri.EO.YJ)  FIJ«FND(A)*!FND!  $J)-FMP!Sin 

ir(YI.MC.YJ)  CALL  S  I  MSON(  him,  ri  I,  SI  ",  S  JP,  .CCI  ,t  ) 

P ETHRN 
ENJ 


FMNCTIIN  Funis) 

COMMON  /ABE/  A,B,C 

FNN-.3 'B'»A228G4*"X'»(-.'*S*S  J*rND(  A»B*(S-E  ») 

RST’IRN 

END 
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1 

Commander 

12 

Defense  Technical  Information 
Center 

ATTN:  DDC-TC 

Cameron  Station,  VA  22314 

US  Army  Materiel  Development 
and  Readiness  Command 

ATTN:  DRCCP 

DRCDE-F 

DRCRE 

2 

Commander 

US  Army  Armament  Research  & 
Development  Cormand 

ATTN :  DRDAR-SEA 

Technical  Library 

Dover,  NJ  07801 

DRCPA-S 

DRCQA 

DRCDE-R 

DRCDE-D 

5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

1 

1 

Commander 

Rock  Island  Arsenal 

ATTN:  Tech  Lib 

Rock  Island,  IL  61299 

Commander 

1 

Comnander 

Harry  Diamond  Laboratories 
ATTN:  DELHD-SAB 

2800  Powder  Mill  Road 

Adel  phi ,  MD  20783 

US  Army  Test  &  Evaluation 
Command 

ATTN:  STEDP-MT-L 

Dugway  Proving  Ground,  UT 

84002 

1 

Commander 

US  Army  Mobility  Equipment 
R&D 

ATTN:  DRDME-0 

Fort  Bel  voir,  VA  22060 

1 

Conmander 

US  Army  Natick  RSD  Command 
ATTN:  DRDNA-0 

Natick,  MA  01760 

2 

Chief 

Defense  Logistics  Studies 
Information  Exchange 

US  Army  Logistics  Management 
Center 

I 

Commander 

US  Army  Concepts  Analysis 
Agency 

ATTN:  DRXMC-D 

Fort  Lee,  VA  23801 

1 

8120  Woodmont  Avenue 

Bethesda,  N)  20014 

Commander 

1 

Reliability  Analysis  Center 
ATTN:  Mr.  I.  L.  Krulac 
Griffiss  AFB,  NY  13441 

US  Army  Aviation  RSD  Command 
ATTN:  DRDAV-BC 

P.0.  Box  209 

St.  Louis,  MO  63166 

1 

Commander 

US  Army  Electronics  RAO 
Command 

ATTN:  DRDEL-AP-OA 

2800  Powder  Mi  1 1  Road 

Adel  phi ,  MD  20783 
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1 

Commander 

US  Army  TRADOC  Systems 

US  Army  Missile  Command 

Analysis  Activity 

ATTN:  DR SMI -C  (R&D) 

ATTN:  ATAA-SL 

Redstone  Arsenal ,  AL  35898 

ATAA-T 

White  Sands  Missile  Range 

1 

Commander 

NM  88002 

US  Arr^y  Troop  Support  & 

Coimtander 

Aviation  Materiel  Readiness 
Command 

US  Army  Tank-Automotive  Command 

ATTN:  DR STS- BA 

ATTN:  DRSTA-TSL 

4300  Goodfellow  Blvd. 

DRSTA-V 

St.  Louis,  MO  63120 

Warren,  MI  48090 

1  Commander 

Air  Force  Armament  Laboratory 
(AFSC)  AFATL/DLMA 
ATTN:  Lt.  James  R.  McClendon 
Eglin  AFB,  FL  32542 


Aberdeen  Proving  Ground 


2  Cdr,  USATECOM 

ATTN:  DRSTE 

DRSTE-CS-A 
Bldg  314 

Dir,  BRL,  Bldg  328 
Dir  BRL 

ATTN:  DRDAR-TSB-S  (ST1NF0  BRANCH) 
Bldg  305 

Dir,  HEL,  Bldg  520 
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